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Abstract 

Since EXOSAT produced the first good quality X-ray light curves of Seyfert galaxies 
there have been several claims of quasi-periodic oscillations (QPOs). None of these 
have withstood repeated analyses and observations. We review some problems con- 
cerning the detection of QPOs. 
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1 Introduction: Why look for QPOs in Active Galaxies? 



Active Galactic Nu clei (AGN) show stro ng, erratic brightness fluctuations 



in the X-ray band ((Lawrence et al.lll987l h At least for the Seyfert galaxies, 



the power spectrum of these variations is a power law-like continuum with a 
steep slope, except in a few cases where the slope can be seen to flatten at 
low frequencies (see e.g. Uttlev et al. 20021 : Markowitz et ah 20( p). Pr ocesses 
with steep power spectra like this are known as red noise ((Press 1978j h 



The X-ray variations in Seyfert galaxies bear more than a p assing resemblance 
to those of X-ray binaries (XRBs; e.g. van der Klisl 19951 ). In addition to a 
relatively featureless red noise continuum, XRBs often show QPOs in the form 
of peaks above the continuum (red noise) power spectrum. By analogy with 
XRBs, AGN may also possess QPOs. For example, assuming a 1/Mbh scaling 
of frequencies we might expect to see analogu es of the high frequency QPOs 
of XRBs at ~ 3 x lO- 3 (M B ff/lO 6 M )- 1 Hz (lAbramowicz et all l2004h . The 
detection of periodic or quasi-periodic variations from a Seyfert galaxy would 
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be a major discovery and could lead to a breakthrough in our understanding 
if the characteristic (peak) frequency could be identified with some physically 
meaningful frequency. 



2 Textbook search methods 



One of th e most popular tools for time series analysis in astronomy is the peri- 
odogram ( Schuster . Il898l ) which is simply the modulus-squared of the discrete 
Fourier transform of the time series (with an appropriate normalisation) . This 
measures the amount of variability power as a function of temporal frequency. 
A concentration of power in a narrow range of frequencies would indicate a 
QPO 1 . The notorious problem with the periodogram of noise processes is 
that the ordinates a re distributed li ke \ 2 (with two degrees of freedom) about 
the true spectrum (|Priestlevl . 1 1 9 8 lh . This means there are large fluctuations 
in the periodogram that can easily be mistaken for genuine peaks (periodic 
signals). Tests for periodic variations against a background of white (flat spec- 
trum) noise tha t take a c count of these intrinsic sa mpling flu c tuatio ns are well 



established: see Fisher (1929) or sectio n 6. 1.4 of |Priestlevl (Il98llh and also 



iLeahv et all (|l983l h lvan der Klisl (|l989h and lPress et al l (|l99fih . If the under 
lying source variability is not white noise the problem is more difficult and 
these methods will yield unreliable results ( van der Klisl 19891 : Benlloch et al. 



2001). 



3 Possible sources of error 



3.1 Systematic/instrumental effects 



With non- imaging detectors (e.g. EXOSAT, Ginga, RXTE) there is always the 
possibility the light curve of a particular AGN is contaminated by other bright 
X-ray sources within the field of view. This possibility was realised in the case 
of NGC 6814 which had been identified as the best candidate for a Seyfert 
galaxy with an X-ray period based o n non- imaging EXOSAT observations 
( Mittaz fc Branduardi-Ravmont 19891 ). Subsequent X-ray imaging observa- 
tions made with ROSAT showed the periodic oscillations belonged to another 
sourc e (V1432 Aquilae, a Cataclysmic Variable) lo cated 37 arcmin from NGC 



6814 ( Madeiski et al. 19931 : Staubert et al. 19941 ) . It is also possible for sys 



tematic, periodic changes in the detector properties to artificially modulate 

1 In the following we shall denote any strictly or quasi-periodic oscillation with the 
generic label "QPO." 
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the observed light curv e . An example is the EXOSAT data of NGC 5548. 
Papadakis fc Lawrence ( 1993) reported significant periodic ities in this mod- 



erate l uminosity Seyfert but a reanalysis of the same data bv lTagliaferri et al. 



(1996) demonstrated that one observation was affected by detector problems. 
The likelihood of these problems is greatly reduced with the current generation 
of low-background, imaging X-ray missions. 



3.2 Inappropriate null hypothesis: Assuming white noise 



Since the days of EXOSAT we have known that the X-ray variations of Seyfert 
galaxies have an intrinsically red power spectrum. This means the periodogram 
is not flat in expectation (as for white noise) and is generally not known a 
priori. As a result i t is meaningless to apply the standar d 'periodogram peak' 
tests f)Fisherlll929[ iLeahv et al.lll983[ | Press et aDll996f ). These all assume a 
null hypothesis of independent (white), Gaussian noise. Thus a high rejection 
probability for the null hypothesis obtained from any of these tests merely 
rejects white noise as the sole origin of the variance. It does not reject the 
more likely possibility that there are intrinsic red noise variations in the light 
curve (i.e. non- white, but also not periodic variations). Even a slight amount 
of intrinsic red noise variance will raise the spectrum above the white noise 
expectation and lead to spurious detections. 



For example, Iwasawa et al.l (|l998f ) discussed a possible periodicity in IRAS 
18325 — 5926 from a long (~ 5 day) ASCA observation. The significance of 
the candidate periodicity was estimated by testing against the hypothesis of 
purely Poisson noise. However, once the red noise variations intrinsic to the 
source were included in the periodogr am fitting the ca ndidate periodicity was 
no longer significant at the 95% level ( Vaughan 2005j ). Other Seyfert galaxies 
for which candidate QPOs were report ed without rigorously testin g an origin 
in red noise include RX .10437.4- 4711 (lHarnern fc Marshall Ill996h . MCG-6- 
30-15 dLee et al. Il2000h. Mrk 766 (|Boller et al.ll200lh and IRAS 13224 - 3809 



(Pfcffcrkorn et al. 200. 



m. 



3.3 Subtle power spectral effects 



There are a number of ways that one can account for red noise in a light curve. 
But this is often not a trivial task and there are many pos sibilities for oyeres- 
timating the significance of candidate QPOs. For example. iFiore et al. 1 (|l989h 
reported candidate QPOs in NGC 4151 based on three EXOSAT observations. 
The authors claimed high (> 99%) significance peaks in the periodograms 
after the data were "pre-whitened" by detrending the light curves using a 
parabola. However, no account was made of the uncertainties involved in this 
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Fig. 1. Periodogram of the 1986/060 EXOSAT observation of NGC 4151 measured 
from the background-subtracted 1 — 9 keV ME data binned to 120 s (obtained 
through the EXOSAT data archive at HEAS ARC / GSFC) . The histogram shows 
the raw data and the solid curve marks the Maximum Likelihood fit using a model 
comprising a power law (red noise) plus constant (Poisson noise). The dashed line 
represents a (minimum) 95% confidence limit for the detection of a QPO calculated 
as a factor j t = — ln[l — (1 — e) 1 /™] above the model, where n is the number of 
frequencies and e is the 'false alarm probability'. This quick calculation slightly 
underestimate s the threshold b ecause it does not account for uncertainties in the 
modelling (see IVaiighan 11200.4) The arrow indicates the position of the candidate 
QPO discussed bv lFiore et al.1 il98flh . 



procedure, which would undoubtedly leave some residual red noise component 
in the spectrum. If the continuum power at low frequencies (where the QPOs 
occurred and where the red noise will be strongest) was underestimated by 
only a small fraction this would drastically innate the significance of the can- 
didate QPOs. Figure 1 shows a re-analysis of one of these datasets. Explicitly 
fitting the red noise and Poisson noise components of the spectrum allows us 
to place more realistic limits on the significances of any QPOs, which now fall 
below the 95% confidence detection threshold. 



The above example is one where the uncertainty in modelling the aperiodic 
power spectrum (i.e. the intrinsic red noise and instrumental Poisson noise) led 
to an artificially enhanced periodicity significance. In general, not accounting 
for the uncertainty in the null hypothesis model (the continuum) often tends 
to overestimate the significance of any peaks in the residual ("whitened") 
periodogram. For example, in their r e-analysis of the EXOSAT observations 
of NGC 5548. lTagliaferri et al. I Jl996) found that in every case the significance 
of the candidate QPO was lower than previously reported (< 95%) once due 
allowance was made for the uncertainties in modelling the spectrum. 

It is also possible to overestimate the significance of a per iodogram peak 
when the periodogram is oversampled. lHalpern et al. ( 2003| ) identified pos- 
sible QPOs in three Seyfert 1 galaxies (Ton S180, RX J0437.4 - 4711, 1H 
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0419-577) based on a uniform analysis of all the available long EUVE observa- 
tions. The significances of the periodiciti es were estimated using a Monte Carlo 
procedure similar to that discussed bv iBenlloch et aL ( 2001 ). However, the 
periodograms were oversampled (computed at frequencies between, as well as 
at, the standard Fourier frequencies). This increases the number of f requencies 



searc hed and therefore also the chance of seeing a spurious signal (|Vaughan 
2005b . 



Finally there are cases where genuinely significant power s pectral features have 
been d etected but were not caused by periodic oscillations. IPapadakis fc Lawrence! 
(1995) reported a broad bump in the power spectrum of NGC 4051 (based on 
EXOSAT observati ons). However, recent observations with R OSAT, RXTE 
and XMM- Newton flOreen et al. Ill999t lM c Hardv et al~ll2004h have provided 
much tighter constraints on the power spectrum and have shown that the un- 
derlying continu um spectrum of NGC 4051 br eaks at a similar frequency to 
that reported by Papadakis fc Lawrence! ( 19951 ) . The broad bump in the spec- 
trum reported from the EXOSAT observation is therefore likely to be due to 
the break in the continuum spectrum rather than an additional QPO feature. 



4 Conclusions 



A major source of 'false alarm' QPOs is the use of an inappropriate contin- 
uum power spectrum. Arbitrarily assuming a white noise spectrum, or red 
noise with an ad hoc slope, will alter the apparent strength (significance) of 
any peaks in the residual periodogram. The continuum spectrum must be 
modelled explicitly. Furthermore, one must account for uncertainties in the 
continuum modelling. Failure t o include these may result in an overly op- 
timistic signific ance calculation ( Israel &: Stella 19961: Tagliaferri et al. 19961: 
Vaughan lEoolh . 



The above arguments apply not only to searching for peaks in the periodogram 
but also for searches based on epoch folding, whereby one bins the time series 
into phase bins at a test period. At the correct period the periodic variations 
will sum while any background noise will cancel out, thus revealing the profile 
of the periodic pulsations. However, the background noise will only cancel out 
if it is temporally independent, i.e. white noise. Again, the presence of any 
underlying red noise variations will produce unreliable res ults unless checked 
with e.g. Monte Carlo simulations 2 (jBenlloch et al.l l2001). 



2 A further word of caution is required here. For the correct application of Monte 
Carlo methods it is crucial to use a random noise generator that correctly produces 
the statistical aspects of red noise and the observational sampl ing, otherwise the 
significances may again be badly biased ( Timmer Koniel 1995h . 
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